Low-frequency incommensurate magnetic response in strongly correlated systems 
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It is shown that in the t-J model of Cu-O planes at low frequencies the dynamic spin structure 
factor is peaked at incommensurate wave vectors (| ± <5, i), (i, | ± S). The incommensurability 
is connected with the momentum dependencies of the magnon frequency and damping near the 
antiferromagnetic wave vector. The behavior of the incommensurate peaks is similar to that observed 
in La2- a: (Ba,Sr) a ;Cu04+ y and YH^Cu^Or^y-. for hole concentrations 0.02 < x < 0.12 we find that 
5 is nearly proportional to x, while for x > 0.12 it tends to saturation. The incommensurability 
disappears with increasing temperature. Generally the incommensurate magnetic response is not 
accompanied by an inhomogeneity of the carrier density. 

PACS numbers: 71.10.Fd, 74.25.Ha 



One of the most interesting features of the inelastic 
neutron scattering in lanthanum cuprates is that for hole 
concentrations x ^ 0.04 the spin fluctuation scattering 
remains commensurate: it is peaked at the antiferromag- 
netic wave vector Q = (5,5) in the reciprocal lattice 
units 2n/a with the lattice period a. However, for larger 
x and low energy transfers the spin fluctuations become 
incommensurate with scattering peaks shifted to the four 
symmetric positions (\ ± 5, ^) and (|, 5 ± S) Q. For 
x < 0.12 the incommensurability parameter S is approx- 
imately equal to x and varies nearly linearly with the 
temperature of the superconducting transition T c ■ For 
larger x the parameter saturates near the value S ~ 0.12. 
The incommensurate response was observed both below 
and above T c Q. There are no clear experimental in- 
dications that in L^-zSr^CuC^ the dynamic magnetic 
incommensurability is connected with lattice distortions. 
However, for x « 0.12 when the static incommensura- 
bility is observed for temperatures T < T c a remarkable 
softening of longitudinal sound waves and soft optical 
phonons were observed Q. This softening is a mani- 
festation of the incipient structural transition from the 
low-tempcraturc orthorhombic to the low-temperature 
tetragonal phase and the appearance of the static incom- 
mensurability is connected with the pinning of the dy- 
namic incommensurate spin correlations by this lattice 
instability |4j. 

The experimental data suggest that the magnetic in- 
commensurability emerges due to doped holes. Sev- 
eral approaches have been used to clarify the underlying 
mechanism. As possible candidates the Fermi surface 
nesting mean- field spiral Q and stripe R| phases, 
and phenomenological marginal Fermi liquid ja| were con- 
sidered. In works based on the Fermi surface nesting 
the difference in the low-frequency magnetic response in 
La2- 2: Sr :r Cu04 and YBa2Cu307_ a was attributed to the 
difference in their Fermi surfaces calculated in the mean- 
field approximation. However, recently the incommen- 
surability which is analogous to that observed in lan- 



thanum cuprates was also detected in YBa2Cu307_ y 
This result indicates that the low-frequency incommen- 
surate magnetic response may be a common property of 
cuprates. This is especially interesting because the over- 
all magnetic response of these two classes of crystals dif- 
fers essentially. As known, in YBa2Cu307_ a and some 
other high T c cuprates the resonance peak is observed 
[rfT | at somewhat larger frequencies u> — 20 — 40 meV 
than the frequencies of incommensurate response, while 
no such peak was detected in lanthanum cuprates 
In works based on the stripe mechanism the magnetic 
incommensurability is connected with the appearance of 
the charge density wave. In neutron scattering such a 
wave is observed only in the low-temperature tetrag- 
onal or the low-temperature less-orthorhombic phases 
where the wave is stabilized by the corrugated pattern 
of the in-plane lattice potential Among lanthanum 
cuprates La2- x Ba 2; Cu04 and La2-j,- x Ndj,Sr 2 .Cu04 have 
these phases. In I^-zSr^CuCU in the low-temperature 
orthorhombic phase no indications of the charge density 
wave were observed in neutron scattering 13] . 

In this letter we demonstrate that the low-frequency in- 
commensurate response can be interpreted based on the 
t-J model of Cu-0 planes. We use the gene ral form of the 
spin spectral function B(kcu) obtained [l4j in Mori's pro- 
jection operator technique. This approach allows us to 
take proper account of strong electron correlations inher- 
ent in cuprates. In calculating the magnon damping the 
known shape of the spin-polaron band |lfj| which forms 
the Fermi surface at low x was used. The mechanism 
of the appearance of the incommensurability is the fol- 
lowing: at low frequencies B(kuj) is reduced to the quo- 
tient of the magnon damping T(k) and the fourth power 
of the magnon frequency lo^. Near Q this frequency 
is an increasing function of |k — Q|. Thus, B(kui) will 
be peaked at incommensurate momenta if r(k) is also 
an increasing function of k— Q|. This takes place in- 
deed because the magnon damping connected with the 
decay into an electron-hole pair is approximately propor- 
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tional to S = J^k'H — n f(£k+k' — M)] n F(£k' — m) where 
ni?(w) = [exp(w/T) + £k is the hole dispersion and 
H the chemical potential. Due to the approximate pe- 
riodicity of £k with the period Q (connected with the 
short-range antiferromagnetic order) S = for k = Q 
and grows with growing |k — Q|. 

The spin spectral function, connected with the re- 
tarded spin Green's function D(ka>) = ((sl\s z ,,)) by the 
relation B(kuj) = — 7r Im.D(ku;), reads [14| 



7r 1 cijr(ka;) 



[uj 2 - f-^Qikiu) - UJ 2 } 2 + /-2 Cl) 2 r2 ( ka; ) ' 



(1) 

where s£ is the z-component of the spin-i operator de- 
scribing localized Cu spins, / = 4JCi(7k — 1) with 
7k = j ^ a exp(ika), a are four vectors connecting near- 
est neighbor sites in the square lattice of the 2D t-J 
model, J is its exchange parameter, C\ = 2(s^s^ +a ) is 
the spin correlation on the neighbor sites n and n+a with 
the statistical averaging denoted by the angular brackets. 
In the normal state the magnon frequency and the damp- 
ing are given by the relations 



r(kw) = 



16aJ 2 |C 1 |(l-7k)(A + l + 7 k), 



327rt 2 J 2 1 - exp(w/T) 



(2) 
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x [1 - iif(ui)]riF(w + uji - a7 2 )-B(k 2 W2) 
x A(kiwi)A(k + ki -k 2 ,uj + uj 1 -uj 2 ), (3) 

where t is the hopping parameter of the t-J model, N is 
the number of sites, a is the vertex parameter |1(t| which 
is of the order of unity, gkk^ = (7k 2 + i)(7k 3 -ki - 
7k x - 7k+ki-k 2 + Tk+kJ, n B (uj) = [exp(u/T) - l]" 1 , 
and A(kui) = — 7r _1 Im ((akerlaj^)} is the hole spectral 
function with the operator creating a hole with the 
spin projection a — ±1. In E q.pt , the parameter A de- 
scribes the spin gap at Q (see [l6| and references therein) . 
The correction 0(kw) to the magnon frequency can be 
obtained from r(kw) and the Kramers-Kronig relation. 
The dynamic structure factor measured in neutron scat- 
tering experiments is connected with B{kui) by the rela- 
tion SQvuj) = tt[1 - exp(-w/T)]- 1 J B(kw) [3- 

If the magnons near Q are not overdamped, the domi- 
nant feature of the spectral function JU is the maximum 
near uj — ujq = AJ{2a\C\ | A) 1 / 2 . This maximum, the 
resonance peak, is observed in YBa2Cu307_j, and some 
other high T c cuprates [ljj. Notice that its frequency 
uj = 20 — 40 meV is close to the frequency of the spin gap 
wq |l4l Il8j . In lanthanum cuprates the magnons near 
Q are apparently overdamped even in the underdoped 
region and due to the small superconducting gap the res- 
onance peak was not observed in the superconducting 




FIG. 1: The Brillouin zone of the square lattice. Solid curves 
are two of four ellipses forming the Fermi surface of the band 
© at small x. Dashed lines are the Fermi surface contours 
shifted by ±(k — Q). Regions of ki and k + ki contributing 
to the damping Q are shaded. 



state either 0, 0] . In this case the frequency depen- 
dence of the spectral function is determined by r(kw) 
rather than the resonance denominator in Eq. l[T]l. In 
lanthanum cuprates the incommensurate response is ob- 
served for frequencies of the order of several millielec- 
tronvolts. For such frequencies Eq. can be essentially 
simplified, 



B(\slj) 



- 1 wr(kw) 



(4) 



As seen from Eq. (J2J, Wk ~ [ujq + c 2 (k — Q) 2 ] 1 / 2 near 
Q and therefore in Eq. (@J u)^ decreases rapidly with 
increasing |k — Q|. In this case for the function B(ku>) 
to be peaked at incommensurate momenta near Q the 
numerator of the fraction Q , the magnon damping, has 
to be a growing function of |k — Q|. 

To make sure that the damping may really have such 
a behavior let us consider the case of low hole concentra- 
tions when in Eq. Jjyl the hole and spin spectral functions 
can be approximated as [ljj 



A(kw) 



{UJ - £ k + M), 



2 V a VA + l+7k 
x[5(uj-uJk) -S(uj + uj k )], 



(5) 



where <j> = jr fl+x) and the dispersion of the spin-polaron 
band [ll El 

e k /f = -2.3893 + 0.055358[cos(2fe x ) +cos(2fe 9 )] 

+0.17857 cos{k x )cos(k y ). (6) 

The approximate formula JfjJ) reproduces the main known 
features of the spin-polaron band [lij . In particular, it 
has minima at (±7r/2, ±7r/2). The bottom of the band, 
— 2.5i, and its width, OAt = 2 J, are set to be close to the 
values obtained in spin-wave and exact diagonalization 
calculations for J = 0.2t. This parameter ratio is close 
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FIG. 2: The dynamic structure factor vs. wave vector along 
the edge (solid line) and the diagonal (dashed line) of the 
Brillouin zone for x = 0.06, T = O.Olt, and lj = 0.004i. 

to that observed in hole-doped cuprates [l^. Also the 
shape of the band © and its Fermi surface are similar 
to those observed in oxychlorides and in L^-zSr^CuCU 
at light and moderate doping 0, |2£j|. For low x the 
Fermi surface of the band © consists of four ellipses 
centered at (±tt/2, ±tt/2). Two of them are shown in 
Fig. n Substituting Eq. JSJ in Eq. @ and taking into 
account that for k « Q and small u> the region around Q 
makes the main contribution in the summation over k2 
we find 

327rf 2 J 2 2 1 - exp(^/T) 

L[kU) PP u 

x X^kkiQl 1 - npie-k, - fi)}n F (e- k+ - kl - fJ,) 

ki 

x 2J B(\l 2 , u + £ki - ^k+kj -k 2 ) 

k 2 

X n B (uJ + £ki - Ek+kj-ka)- (7) 

The dependence T(k) is mainly determined by the sum 
s = E kl SkkiQl 1 ~ n F (e kl - ^)]n F (e k+kl - fi). For low 
temperatures the regions of ki and k + ki contributing 
to the sum are shaded in Fig. The damping is roughly 
proportional to the area of these regions. Thus, for low 
T and x the damping grows with growing |k — Q| which 
shifts the peaks in B(kuj) to incommensurate positions. 

Key elements for the validity of the above conclusion 
are the spin spectral function in Eq. which is strongly 
peaked near Q and an approximate periodicity of the 
hole dispersion with the period Q. Due to the doubling of 
the elementary cell in antifcrromagnctically ordered crys- 
tals both these conditions have to be fulfilled in lightly 
and moderately doped cuprates which are characterized 
by the short-range antiferromagnetic ordering. A small 
damping of hole states is also necessary (the incommen- 
surability was suppressed in the self-consistent calcula- 
tions | lq due to a sizeable hole damping introduced for 
stabilizing the iteration procedure). 

We used Eqs. © _ © an d parameters a, Ci, and A ob- 
tained self-consistently for calculating the dynamic 
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FIG. 3: (a) The dynamic structure factor vs. wave vec- 
tor along the edge for x = 0.015 (solid line), 0.06 (dashed 
line), and 0.12 (dotted line). For all curves T = O.Olt and 
lo = 0.004f. Inset: the structure factor measured 0, l2l| in 
La2-iSr^Cu04 for x = 0.03, T = 2 K, lu = (filled squares 
and solid line), x = 0.06, T = 25 K, uj — 2 meV (triangles 
and dashed line), and x = 0.12, T = 31 K, uj = 2 meV (open 
squares and dotted line). Fitted curves in the inset are from 
0, l2l| . (b) The incommensurability parameter 8 vs. x for 
T = O.Olt and uj — 0.004£. Inset: experimental data 0] for 
La2- a: Sr a: Cu04. Connecting lines are a guide to the eye. 



structure factor S(kuj). We considered lattices up to 
120x120 sites with the parameters J — 0.2t. Figure |3 
demonstrates the dynamic structure factor calculated 
along edge and diagonal of the Brillouin zone near Q. 
The parameters chosen are close to those used in exper- 
iments: for t = 0.5 eV, T w 58 K and oj = 2 meV. In 
agreement with experimental observations Q, 0, Q peaks 
along the edge are more intensive than those along the 
diagonal. 

Figure [21 shows the concentration dependence of the 
peak position along the edge. Starting from x ks 0.02 
there appears the incommensurability - instead of one 
peak at Q four peaks are observed at (|, | ± S) and 
(i ±6, g). For x < 0.12 the incommensurability param- 
eter S grows nearly linearly with growing x and tends to 
saturation for x > 0.12. As seen from the insets in Fig.|3 
this behavior is similar to the experimental observations 
except that in experiment the incommensurability arises 
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FIG. 4: The dynamic structure factor vs. wave vector along 
the edge for x = 0.04, u = 0.004t and T = O.Oli (solid line), 
0.02t (dashed line), and 0.04t (dotted line). Inset: the struc- 
ture factor in Lai.g3Sro.o7Cu04 measured |22| for u) — 4 meV, 
T = 23 K (filled squares and solid line), T = 100 K (trian- 
gles and dashed line), and for uj = 8 meV, T = 293 K (open 
squares and dotted line). Fitted curves in the inset are from 



at x k, 0.04 and the parameter S is approximately 1.5 
times larger. In our calculations its concentration de- 
pendence is mainly determined by the variation of the 
spin-gap parameter A in Wk in Eq. @. For x < 0.12 this 
parameter grows with x and saturates at x w 0.12 [lfij ]. 

The temperature variation of the peaks is shown in 
Fig. 0] With growing temperature S decreases and finally 
the incommensurability disappears. As seen in the inset 
of Fig. the analogous behavior of the dynamic struc- 
ture factor is observed experimentally [2l|22|. Notice that 
using Eq. J^J we neglected completely the damping in the 
hole spectral function. As mentioned, this damping de- 
creases or even suppresses the incommensurability. The 
thermal growth of the damping hastens the decrease of 
the incommensurability with temperature. 

It should be emphasized that at least in the case when 
magnons near Q are not overdamped the low-frequency 
incommensurability in B(huj) influences only weakly the 
hole spectral functions. The solution of the self-energy 
equations .1(5] for the t-J model indicates that in these 
conditions the density of holes remains homogeneous 
as it is apparently observed in the low-temperature or- 
thorhombic phase in L^-xSr^CuC^ 13]. 

In summary, it was shown that for hole concentrations 
x > 0.02 the low-frequency magnetic response of the two- 
dimensional t-J model is incommensurate. The incom- 
mensurability is the consequence of the local minimum 
in the momentum dependence of the magnon damping. 
The minimum is located at the antiferromagnetic wave 
vector Q = (5,5) in the reciprocal lattice units and is 
connected with the growth of the phase space for the 
magnon damping into the electron-hole pair with dis- 



tance of the wave vector k from Q. For such behavior 
the spin spectral function has to be strongly peaked near 
Q and the hole dispersion has to be approximately pe- 
riodic with the period Q. Both these conditions are ful- 
filled in crystals with short-range antiferromagnetic order 
to which lightly and moderately doped cuprates belong. 
In agreement with experiment for x < 0.12 the incom- 
mensurability grows nearly proportional to x and tends 
to saturation for x > 0.12. This dependence is connected 
with the concentration dependence of the spin gap at 
Q. Also in agreement with experiment the incommen- 
surability decreases with increasing temperature. In the 
considered mechanism the magnetic incommensurability 
is not accompanied by the inhomogeneity of the carrier 
density. 
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